Centrality dependence of charged hadron transverse momentum spectra in d+Au 

collisions at y/s NN = 200 GeV 
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We have measured transverse momentum distributions of charged hadrons produced in d+Au 
collisions at ^/s NN = 200 GeV. The spectra were obtained for transverse momenta 0.25 < pr < 
6.0 GeV/c, in a pseudorapidity range of 0.2 < r\ < 1.4 in the deuteron direction. The evolution 
of the spectra with collision centrality is presented in comparison to p + p collisions at the same 
collision energy. With increasing centrality, the yield at high transverse momenta increases more 
rapidly than the overall particle density, leading to a strong modification of the spectral shape. This 
change in spectral shape is qualitatively different from observations in Au+Au collisions at the same 
energy. The results provide important information for discriminating between different models for 
the suppression of high-pT hadrons observed in Au+Au collisions. 

PACS numbers: 25.75.-q,25.75.Dw,25.75.Gz 



The yield of charged hadrons produced in collisions 
of deuterons with gold nuclei at an energy of y/s NN = 
200 GeV has been measured with the PHOBOS de- 
tector at the Relativistic Heavy Ion Collider (RHIC) 
at Brookhaven National Laboratory. The data are 
presented as a function of transverse momentum (pt) 
and collision centrality. The goal of these measurements 
is to study the modification of particle production due to 
initial state effects in the nuclear medium, in comparison 
to nucleon-nucleon collisions at the same energy. Mea- 
surements from proton-nucleus (p+A) reactions at lower 
collision energies have found that the cross-section for 
hadron production at px of 1.5 to 5 GeV/c rises faster 
than the nuclear size A p|. This observation, which is 
commonly called the Cronin effect, has been described as 
the result of initial state multiple scattering, leading to 
a broadening of the pr distribution . 

The present interest in repeating these measurements 



at higher energies is motivated by results from Au+Au 
collisions at ^/s NN — 130 and 200 GeV. In these collisions, 
the expected scaling of hadron production with the 
number of binary nucleon-nucleon collisions at px of 
2-10 GeV/c is strongly violated 0, 0, H 0. This 
effect had been predicted as a consequence of the energy 
loss of high-pT partons in the hot and dense medium 
formed in Au+Au collisions 0. The interpretation of 
the Au+Au data relics on the understanding of initial 
state effects, including gluon saturation J5], which can 
be investigated with the d+Au data presented here 0. 
Similar measurements are reported in [Tol ITlj . By 
studying the spectra as a function of collision centrality, 
we can control the effective thickness of nuclear matter 
traversed by the incoming partons. 

The data were collected using the PHOBOS two-arm 
magnetic spectrometer |T^ . The spectrometer arms are 
each equipped with 16 layers of silicon sensors, providing 
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charged particle tracking both outside and inside the 
2 T field of the PHOBOS magnet. Additional silicon 
detectors used in this analysis are the central single- 
layer Octagon barrel detector and the three single-layer 
forward Ring detectors located on either side of the 
interaction point. 

The primary event trigger (Level 0) was provided by 
two sets of 16 scintillator counters ("Paddle counters") 
covering pseudorapidities 3 < \r}\ < 4.5. In addition, 
two higher level trigger conditions were used. Collisions 
close to the nominal vertex position z v t x — along the 
longitudinal (z) direction were selected using the time 
difference between signals in two rings of ten Cercnkov 
counters. The counters covered —4.9 < ?/ < —4.4 and 
3.6 < r) < 4.1, respectively. For part of the data set, 
a further online selection of events was accomplished 
using two arrays of horizontally segmented scintillator 
hodoscopes. One array was positioned immediately 
behind the spectrometer and the other at a distance of 
~ 5 m from the interaction point. In combination with 
the known vertex position, a spatial coincidence from 
the hodoscopes was used to trigger on high-py particles 
which traverse the spectrometer. 

Previous measurements of p+A collisions have shown 
a strong dependence of the multiplicities and momentum 
distributions of produced particles on the size of the tar- 
get nucleus For d+Au collisions, this implies the 
importance of characterizing the collision centrality. For 
a given event selection, the centrality can be quantified 
by a number of variables such as the average number 
of participating nucleons (N part ) or the average number 
of binary nucleon-nucleon collisions, (N co u). Neither of 
these variables can be directly measured. Estimates of 
(N par t} or (N co u) for a given event selection are typically 
obtained by matching the distribution of an experimental 
observable, such as the multiplicity in a certain 77-region, 
with the results of a Glauber model calculation and 
detector simulation [bH ]. 



Centrality Selection 


{Np ar t} 


(Ncou) 


Efficiency 


0-20% 


15.5 ± 1.0 


14.6 ±0.9 


82% 


20-40% 


10.9 ±0.9 


9.7 ±0.8 


73% 


40-70% 


6.7 ±0.9 


5.4 ±0.8 


49% 


70-100% 


3.3 ±0.7 


2.2 ±0.6 


14% 



TABLE I: Estimated values for {N par t) and (N co a) and the 
respective systematic uncertainties for the four centrality bins 
used in this analysis. Also shown is the average event selection 
efficiency for each of the bins. The centrality bins are based 
on the signal of multiplicity counters covering 3.0 < \rj\ < 5.4. 



The centrality cuts for this analysis were based on 
the signal (Em ng ) of the three Ring detectors in the 
region of 3.0 < |?7| < 5.4, which is proportional to the 



number of charged particles hitting these counters. In 
the MC simulations, centrality cuts on Em ng were found 
to introduce a bias on the yield in the spectrometer 
acceptance of less than 5% for (N part ) > 3, compared to 
cutting directly on N part . If the N part determination is 
based on a multiplicity region closer to the spectrometer 
acceptance, MC studies, as well as our data, showed a 
significantly larger bias at low and high centralities. For 
each of the four bins in Em ngi (N part ) and (N co u) were 
obtained from a Glauber model calculation using HIJING 
|15| and a full detector simulation which included the 
experimental trigger and event selection efficiency. In 
the HIJING calculations, the default value for the in- 
elastic nucleon-nucleon cross-section of 41 mb was used, 
consistent with previous calculations [l4| . 

The resulting estimates for (N part ) and (N co u) in the 
four Efii n g -based centrality bins are shown in Table I. 
The percentage numbers refer to the fractional cross- 
section in the unbiased HIJING distribution. The 
determination of (N part ) and (N co a) takes into account 
the bias introduced in the centrality by the online and 
offline event selection in the relatively low-multiplicity 
d±Au events. The largest contribution to this bias comes 
from the online vertex trigger, leading to the average 
event selection efficiencies for the individual centrality 
bins shown in Table I. 
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FIG. 1: Invariant yield of charged hadrons, h + h , as a 
function of pr for four centrality bins. Only statistical errors 
are shown. 

Details of the track reconstruction algorithm can be 
found in [|| llG( . To optimize the momentum resolution 
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FIG. 2: The relative yield of d+Au to UA1 p+p data at p T = 
1,2,3 and 4 GeV/c is shown as a function of the integrated 
yield dN/drj for the four Em-ng centrality bins. The ratio 
d+Au/p+p has been normalized to unity at pr = 0.5 GeV/c. 
The triangles indicate the values for p + p. The brackets 
indicate the systematic errors on the relative yield (90% C.L.). 
The systematic error on dN/drj is 12%. 



and minimize systematic errors in the track selection, 
only particles traversing the full spectrometer arms were 
included in the analysis. These particles leave a minimum 
of 12 hits in the silicon detectors. This selection limits 
the usable vertex range to —15 cm < z v ± x < +10 cm. 
Due to the low multiplicity in d+Au collisions, a new 
algorithm for the offline determination of the collision 
vertex was developed, using hit position and energy 
information in the Octagon detector. MC studies show 
a resolution in the beam direction of a v tx z =1.4 cm for 
the most peripheral events and a v t Xz — 0.8 cm for the 
most central events. The transverse position of the event 
vertex was centered at the known position of the beam 
orbit. Unlike in the Au+Au track finding, the vertex 
position information was not included in the initial track 
seed. 

To obtain the invariant yield of charged hadrons, we 
accumulated equal amounts of data with both magnet 
polarities. The transverse momentum distributions for 
each centrality bin were corrected for the geometrical 
acceptance of the detector, the efficiency of the tracking 
algorithm and the distortion due to binning and mo- 
mentum resolution. The procedure for obtaining the 
correction factors was described in 0. The largest 
contributions to the systematic uncertainty come from 
the overall tracking efficiency (5-10% uncertainty) and 
the reduction in overall acceptance due to malfunctioning 
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FIG. 3: Nuclear modification factor RdAu as a function of 
Pt for four bins of centrality. For the most central bin, the 
spectral shape for central Au+Au data relative to p + p is 
shown for comparison. The shaded area shows the uncertainty 
in RdAu due to the systematic uncertainty in (N C oll) and 
the UA1 scale error (90% C.L.). The brackets show the 
systematic uncertainty of the d+Au spectra measurement 
(90% C.L.). 



channels in the silicon detectors (5% uncertainty). The 
corresponding corrections are centrality independent. 
The next largest correction is the px and centrality 
dependent momentum resolution and binning correction. 
The contamination by secondary particles and feeddown 
particles is small, due to the proximity of the tracking 
detectors to the collision vertex and the requirement for 
the reconstructed track to point back to the beam orbit 
to within 0.4 cm. 

In Fig. ^ we present the invariant yield of charged 
hadrons as a function of transverse momentum, obtained 
by averaging the yields of positive and negative hadrons. 
Data are shown for four Em ng centrality bins. The plot 
shows the evolution of overall yield and spectral shape 
with increasing collision centrality. 

The centrality evolution of the spectra can be studied 
in detail in Fig. [21 where we compare our d+Au data to 
results from UAl for p + p collisions at the same energy 
[Tt| . To account for the difference in acceptance between 
UAl < 2.5) and PHOBOS, a correction function was 
determined using PYTHIA 18]. The quantity shown 
on the vertical axis in Fig. [21 is a direct measure of the 
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FIG. 4: Nuclear modification factor -RdAu as a function of 
centrality in four bins of transverse momentum. The brackets 
indicate the point-to-point systematic error, dominated by the 
uncertainty in the number of collisions for each centrality bin. 
The grey band shows the overall scale uncertainty at each pr- 
Systematic errors are at 90% C.L. 



modification of the spectral shape relative to p + p for 
each centrality bin in d+Au and is defined as follows: 
the fitted dN/dpr distribution for each centrality bin 
is divided by the corrected fit to the UAl p + p data 
|l9j | . The resulting distribution is normalized to unity at 
Pt = 0.5 GeV/c. Then the value of the normalized ratio 
at four values of pr from 1 to 4 GeV/c is plotted against 
dN/drj determined by integrating the hadron spectrum 
for each centrality bin. If there were no modification of 
the spectral shape relative to p + p collisions, the ratio 
would be flat at unity. This presentation of the data 
allows an investigation of the evolution of the spectral 
shape with increasing centrality, while eliminating the 
uncertainty associated with the determination of (N co u). 
We observe that the relative yield for all pr regions 
grows smoothly as a function of dN/dr], with the biggest 
increase relative to p + p observed at the largest px ■ For 
all Pt regions, the data extrapolate to the corrected p+p 
fit. 



In Fig. we present the nuclear modification factor 
RdAu as a function of px for each centrality bin, defined 
as 
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dAb 



PP 



d 2 N dAu /dp T dr] 



(Ncou) d 2 a{VAl)pp/dp T dr] 



(1) 



Consistent with our Glauber calculations, we used 
<j™ el = 41 mb. A value of RdAu = 1 corresponds to 
scaling of the yield as an incoherent superposition of 
nucleon-nucleon collisions. For all centrality bins, we 
observe a rapid rise of R d Au from low pr, leveling off 
at pr of w 2 GeV/c. For comparison, we also plot the 
results from central Au+Au collisions at the same energy 
in the lower right panel of Fig-El The average number 
of collisions undergone by each participating nucleon in 
the central Au+Au collision is close to 6, similar to that 
of each nucleon from the deuteron in a central d+Au 
collision. For central Au+Au collisions, the ratio of the 
spectra to p + p rises rapidly up to pr ~ 2 GeV/c, but 
falls far short of collision scaling at larger px , in striking 
contrast to the behavior for central d+Au collisions. 

Predictions for the evolution of R d Au from semi- 
peripheral collisions with (N co u) w 6 to central colli- 
sions were made in two qualitatively different models. 
Perturbative QCD calculations ,20] predict an increase 
in the maximum value of RdAu at pr ~ 3.5 GeV/c by 
15%. In contrast, a decrease in R d Au by 25-30% over the 
same centrality range is predicted in a parton saturation 
model @]. The centrality evolution of RdAu is shown 
in Fig. where the points were obtained from a fit to 
the pt dependence of RdAu in each centrality bin. Our 
data disfavor the prediction from the parton saturation 
model. This suggests that the observed suppression of 
high pt hadrons in Au+Au collisions |3|, Ij, |5|, |(| cannot 
be accounted for by initial state effects that should also 
be present in d+Au collisions. 
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